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We present the first numerical study of the Falicov-Kimball model extended by an on-site hy- 
bridization away from the particle-hole symmetry point. Stable polaronic distortions of the charge- 
density wave (CDW) phase are observed when doping with a single hole. For moderate hole-doping, 
we find phase separation between a hole-rich homogeneous state and the usual CDW state. Ex- 
citonic order is enhanced around the hole-rich region, locally manifesting the global competition 
between this and the CDW phase. Associated with the hole-doping is a fractionalization of the 
electron density between the itinerant and localized electrons. The calculated local density of states 
at the polaron centre reveals this fractionalization to be the result of charge fluctuations induced by 
the hybridization. We propose a scanning tunneling microscope experiment to test our results. 

PACS numbers: 71.28.+d, 71.45.Lr 



The Falicov-Kimball model (FKM) has found use far 
beyond its original role as a minimal model for valence 
transitions 0|. The FKM describes itinerant d-electrons 
interacting via a repulsive contact potential U with lo- 
calized /-electrons of energy ef. this simple physical pic- 
ture is a useful starting point for more complex models of 
mixed- valence (MV) phenomenon Q , heavy fermions , 
electronic ferroelectricity 0, IE Hi an d unconventional su- 
perconductivity Q. The "bare" form of the FKM has 
also attracted much interest: it exhibits a charge density 
wave (CDW) instability at half-filling in any dimension 
D and is considered a basic model of binary alloys @ . 

Apart from the extreme limits D = 1 9j and D — > 
oo [jjj], very little is known about the FKM away from 
half-filling. The numerical results of Lcmariski et al. sug- 
gest a rich 2D phase diagram, with phase separated coex- 
istence of many different orderings ; it would clearly 
be very interesting to assess the stability of these states 
within the context of a more realistic model. The ef- 
fect of dilute hole doping on the CDW is also of inter- 
est: in related models displaying density-wave instabil- 
ities (e.g. the Hubbard model) hole-doping away from 
perfect nesting produces polarons |l2j . Yin et al. have 
studied polarons in the strong-coupling FKM with a weak 
/-electron hopping term, obtaining non-trivial modifica- 
tions of the usual t-J model results @. Of greater rele- 
vance, Liu and Ho have proposed that polarons form in 
a FKM with hybridization, as a precursor to a valence 
transition between states with integral and fractional oc- 
cupation of the /-orbitals Q . 

In this Letter we extend our previous study of the 
FKM with d-f hybridization away from the particle- 
hole symmetry point 01 ■ We perform a Hartree-Fock 
(HF) decoupling of the interaction and exactly diagonal- 
ize the resulting real-space Bogoliubov-de Gennes (BdG) 
eigenequations. We work exclusively in the limit of zero 
temperature where the HF approach is known to be most 
accurate. Doping with a single hole produces a local pola- 



ronic distortion of the CDW phase. Unique to our model, 
we find a fractional partitioning of the electron density 
between the two orbitals localized at the doped hole; this 
is revealed by the local density of states (LDOS) spec- 
tra at the polaron centre. We discuss the origin of the 
fractionalization and a proposed experimental approach 
to observing this novel electronic state. For a moder- 
ate hole-doping, we find that phase separation between 
a hole-rich homogeneous phase and the usual half-filled 
CDW phase is the ground state. The global competition 
between these phases is manifested as a local enhance- 
ment of excitonic ordering. 

We consider an extended FKM for spinlcss fermions: 

i j 

In addition to the usual FKM terms, we include a hop- 
ping t{j for the /-electrons and an on-site hybridization V 
between the d- and /-orbitals 0|. Here we concentrate 
on the case V ^ and t{j = 0; we consider also a finite 
tjj with V — for comparison with previous work 0, [|| . 
The hopping integrals tij and t{j are, respectively, t and 

for nearest neighbours, vanishing otherwise. 

We perform the standard HF decomposition Q of the 
Coulomb repulsion term: TijUj = (nj)r^ + (n^)n^ — 

Ajdjfj — Ajfjdj. Here Aj = (fjdj) is the excitonic aver- 
age at site j. In the absence of a hybridization potential 
V, Aj ^ indicates the excitonic insulator phase [TEI ]: 
Portengen et al. interpreted such a "spontaneous" exci- 
tonic average as evidence of electronic ferroelectricity Q . 
The resulting HF Hamiltonian is diagonalized by the 
canonical transform, 7„ = J2j{u r jdj + vj/j). The quasi- 
particle wavefunction amplitudes, u" and u", are derived 



2 



by solving the associated BdG eigencquations: 

S(#8)(?M?)- < 2 > 

where the components of the Hamiltonian matrix are 
defined as Tiff = -t %] + U{nfoi j} H{/ = -t£. + 

(e f + U(nj))Sij and = (V - UA^. In terms 
of the diagonal basis the order parameters at site j 
are given by (r$ - En=l K"l 2 ; <»»/> = E^i Kf 

and Aj = J2n=l v T u j- We also define the d " and /■ 
electron CDW order parameters 8j = (— l) J ((n^) — |) 

and 6? = - §). Note that in the CDW phase 

sgn(<5^) = — sgn(^J). Our calculations are performed in 
the canonical ensemble at zero temperature and thus the 
sum over n extends over the first N to t occupied quasipar- 
ticle states. 

We solve the BdG eigenequations self-consistently us- 
ing a numerical iteration scheme. Commencing with an 
initializing set of order parameters, we exactly diagonal- 
izc Eq. (j2J) and hence compute new order parameters us- 
ing the obtained quasiparticle wavefunctions. These val- 
ues are then used as an input for the next iteration; this 
procedure is repeated until a desired accuracy is reached. 
The calculations are performed for a N = 24 x 24 lat- 
tice with periodic boundary conditions. To calculate the 
quasiparticle LDOS we use the converged 24 x 24 lattice 
as a superccll in a 10 x 10 array. The d-electron hopping 
integral t defines our energy scale. The results presented 
below correspond to e/ =0 and U = 2.0; for the case of 
a single hole wc take V = 0.1 (t f = -0.1), while V = 0.2 
(t f = -0.2) for the case of moderate hole-doping. 

Limit of a single hole. We dope a single hole into 
the CDW state by computing the order parameters for 
Ntot = N —1 (note N to t = N corresponds to half filling) ; 
the converged results are presented in Fig. ^ The hole 
appears as a localized vacancy in the occupied /-electron 
sub-lattice [see Fig.^a)]. This defect in the periodic HF 
potential U (nj) experienced by the ci-electrons produces 
a distortion of the (i-electron CDW state which extends 
over several lattice constants [Fig. IHb,c)], defining the 
spatial extent of a polaron. Near the polaron, the cx- 
citonic average Aj [Fig. ^d)] shows a Friedel-like oscil- 
lation with its tails extending along the diagonals. The 
enhancement of |Aj| about the /-hole indicates that the 
polaron in the FKM is a local manifestation of the com- 
petition between the global excitonic ordering and CDW 
states [l3j . 

Away from the polaron 6j is uniform to within the 
accuracy of the convergence, indicating the expected sin- 
gle d-electron per unit cell in the CDW bulk. Overall, 
however, we find a total (i-electron population N c i = 
288.1084, as against the case of Ndo = 288 d-electrons 
in the usual half-filled CDW state. A "fraction" AN d = 
Nd — Ndo = 0.1084 of a d-electron is associated with 
the polaron. Correspondingly, the total /-electron pop- 
ulation Nf = 286.8916 such that AN f = N f - N f0 = 




FIG. 1: (color online) The variation of the /-electron den- 
sity (a), the d-electron density (b), the d-electron CDW order 
parameter (c), and the excitonic order parameter (d) in the 
vicinity of a single hole. 



— 1.1084, showing the transfer of charge from the /- to the 
d-orbitals. Although the formation of a polaron about 
the /-hole is indicated by Liu and Ho's work 0], their 
method did not allow them to observe the fractional- 
ization. Our result also differs from the homogeneous 
distribution of fractional charge in the traditional MV 
state 0: we find the fractional charge localized at the 
polaron. 

This fractional charge arises from the mixing of the 
orbital wavefunctions by the hybridization. As a result 
of this mixing, the number of d- and /-electrons are not 
separately conserved. (Note that the total number of 
electrons, N = Nd + Nf, remains constant.) The hy- 
bridization allows electrons in the d-orbitals to tunnel 
into the /-orbitals and vice versa: this implies charge 
fluctuations between the two orbitals leading to the ob- 
served fractional populations. Although true also in the 
half-filled CDW state, it is only with the introduction 
of the inhomogeneities by hole-doping (i.e. the polarons) 
that the fractional charge can be directly observed. Frac- 
tionally charged polarons are therefore unique to multi- 
orbital models such as the FKM considered here; there 
is no analogous effect in single-orbital systems. 

We find that the fractional d-electron ANd (or /- 
electron ANf) is strongly dependent upon U and V: 
ANd increases with both these parameters, although it 
decreases as V is raised above the critical value that 
destabilizes the global CDW phase. According to our 
HF decomposition the excitonic average enters into the 
equations as an effective hybridization potential —UAj. 
In the CDW state, however, this is only non-zero in the 
presence of a finite on-site hybridization [T^ |. As the 
magnitude of Aj is directly proportional to both U and 
V, the observed behaviour of ANd is easily understood. 
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FIG. 2: (color online) The d- and /-electron LDOS spectra 
at the polaron centre and in the bulk for (a) V = 0.1 and (b) 
— —0.1 (V = 0). In both panels the bulk d- and /-spectra 
are the blue and green lines, respectively; the polaron d- and 
/-spectra are given by the black and red lines, respectively. 




The fractionalization can be explicitly visualized by 
computing the LDOS both at the polaron site and in 
the CDW bulk. From our converged solutions for the 
order parameters, we calculate the LDOS of the d- and 
/-electrons at site j: 



FIG. 3: (color online) The variation of the d-electron DOS 
p d (co) across the lattice for (a) lo = —0.131 and (b) lo — 
0.168. In (c) and (d) we plot the /-electron DOS p f (lo) for 
the energies corresponding to (a) and (b). 



p%u) = Y,K?5<*{u n -u), (3) 

n 

= £KT^K-^). (4) 

n 

Here 8 a {ui) is the Lorentzian with phcnomcnological life- 
time broadening characterized by the parameter a. In 
our plots of the LDOS spectra we always take a = 0.05. 
We plot the LDOS spectra in Fig. Eta); in the accompa- 
nying Fig. Elb) we plot the same spectra for the FKM 
extended by a finite /-hopping integral = —0.1. 

The effect of the hybridization on the d- and /- 
wavefunctions can be clearly appreciated in Fig. EI a )- 
The off-diagonal terms TL^ in Eq. produce quasipar- 
ticles with mixed nature, as evidenced by the coincidence 
of the peaks in the d- and /-electron spectra. The sit- 
uation is qualitatively different for the FKM extended 
by a finite /-hopping integral [Fig. Elk)] . Essentially a 
two-band Hubbard model, in this case iV^ and Nf sepa- 
rately commute with the Hamiltonian (i.e. Ttf^ = 0) and 
the quasiparticles retain purely d- or /-character. This 
can be seen in the polaron LDOS spectra as two distinct 
peaks for the d- and /-electrons within the CDW gap. 

In spite of the orbital-mixing, the LDOS spectra of the 
two species is still an experimentally accessible quantity; 
for example, we can imagine a scanning tunneling mi- 
croscope (STM) experiment where the tunneling matrix 
only couples the tip to the d-orbitals. The polaron could 
then be identified by examining the spatial distribution 
of the d-electron LDOS at the energies corresponding to 
the two intra-gap peaks in Fig. EI a ). This will distin- 
guish between the fractionally charged polaron and the 
usual Hubbard-like polaron: in our system a peak in the 
d-electron LDOS will be observed at both energies, as 



opposed to only a single peak at the higher energy when 
V = Q. 

Using Eq. © we calculate the variation of pj(ui) across 
the lattice for ui = —0.131 and tu = 0.168. The results 
are presented in Fig. O^a) and (b), respectively; the cor- 
responding LDOS for the /-electrons are given in (c) and 
(d). As expected, for both species we find distinct peaks 
in the LDOS at the polaron site for the two intra-gap 
energies. The d-electron LDOS Fig. E^a) for the (pri- 
marily /-character) lj = —0.131 peak displays the diag- 
onal star-shape seen before in the plot of Aj [Fig. Hfd)]. 
Since the enhancement of | Aj | near the polaron indi- 
cates strong d-f mixing, this structure in the LDOS is 
a defining signature of the fractionally charged polaron. 
A much weaker diagonal structure is present in the /- 
electron LDOS Fig. Etc). The same effect occurs in the 
lo = 0.168 d-electron plot [Fig. EJb)], but it is somewhat 
obscured by the higher d-electron peak; no such structure 
is visible in the associated /-electron image. 

Finite hole doping. To study the effects of moder- 
ate hole-doping, we solve the BdG eigencquations for 
Ntot — N — 10. The results for V = suggest a large va- 
riety of possible charge-orderings [ilT | ; to find the ground 
state we therefore consider a number of different initial 
order parameter sets. These sets range from vertical 
and diagonal stripe orderings to a random distribution of 
quasiparticle weight across the lattice. Since we work in 
the zero temperature limit, we test the relative stability 
of the converged results by calculating the total energy 
E = En E n~ UY,j{{nj)(n*) - l A jf)- We present the 
variations of Sj for a representative sample of these ini- 
tializing sets in Fig.^a-d); these figures arc labelled with 
decreasing energy from E/N = —1.5554 for Fig. Efa) to 
E/N = -1.5586 for Fig.^d). Fig. He) gives the varia- 
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FIG. 4: (color online) Variation of 8f for different initializ- 
ing order parameter sets: (a) period- 12 horizontal stripe; (b) 
CDW; (c) period-24 diagonal stripe; (d) random, (e) Vari- 
ation of Aj corresponding to (d). (f) Variation of Sf for 
t s = -0.2 (V = 0) using same initializing sets as for (b). 

tion of Aj for the lowest-energy solution. In Fig.|3Jf) we 
present Sj for finite t* . 

Our results clearly demonstrate the tendency for the 
/-holes to cluster. Such phase separation is well known 
in the V = FKM, bein g re ported in the 2D model at 
all coupling strengths (l"ltll6| . Our observation of /-hole 
clustering again confirms the relevance of the V = re- 
sults as a guide to the physics of the extended model 01 • 
An attractive interaction between /-holes in a closely- 
related FKM has recently been cited as a possible mech- 



anism for the unconventional "valence fluctuation super- 
conductivity" observed in several Ce compounds Q ■ Al- 
though our results cannot confirm such a mechanism, 
they lend support to this scenario. With the inclusion 
of a finite /-hopping, we see from Fig. 0Jf) that there 
appears to be some tendency for the phase separation to 
transform into a diagonal stripe phase. This tendency 
can be understood by the observation that the spinlcss 
FKM model with v = t can be mapped onto the Hub- 
bard model, where the stripe phase is always obtained 
upon finite hole doping p*7j . 

The fractionalization effect is present in Fig.0Ja-d); in 
particular, for Fig. g^d) we find AN d = 4.0038, indicat- 
ing the transfer of four electrons from the /-orbitals. As 
before, we find no evidence of this transfer in Fig. 0Jf) 
where = Ndo = 288. Also anticipated by our single- 
hole results, we find enhancement of the excitonic order 
parameter Aj|, mainly at the boundary of the /-hole 
cluster [Fig. 0(e)]. This again locally evidences the com- 
petition between the excitonic and CDW orders. 

In conclusion, we have investigated the FKM extended 
by the d-f hybridization away from the half-filling sym- 
metry point. Throughout this report, we have compared 
and contrasted our results with the extended FKM pro- 
posed in We have studied the formation of polarons 
when a single hole is doped into the CDW state. As- 
sociated with the polaron, we have discovered a frac- 
tional partitioning of the electron density between the 
d- and /-orbitals. Such an effect is unique to our ex- 
tended FKM: analysis of the LDOS spectra at the po- 
laron centre reveals the origins of the fractionalization as 
the hybridization term. We propose a STM experiment 
to verify our results. For moderate hole-doping, we have 
demonstrated that phase separation into hole-rich homo- 
geneous and CDW states is energetically favorable. In 
both doping limits, enhancement of the excitonic order- 
ing is found around the hole- rich region, locally manifest- 
ing the global competition between the two phases 0] . 
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